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Abstract 

i?-parity violating supersymmetry in a Minimal Flavor Violation paradigm can produce 
same-sign dilepton signals via direct sbottom-LSP pair production. Such signals arise when 
the sbottom hadronizes and the resulting mesino oscillates into an anti-mesino. The first 
bounds on the sbottom mass are placed in this scenario using current LHC results. 

I. INTRODUCTION 

The 2011 and 2012 data from the Large Hadron Collider (LHC) place severe constraints 
on natural R-parity conserving models of supersymmetry (SUSY) [TJ-[3]. While such models 
are not excluded by the data, if they are to solve the hierarchy problem of the Standard 
Model (SM), they are forced to have either non-generic spectra where only third-generation 
squarks are light [1HH] or nearly degenerate particles, either in the form of stealth SUSY [9] 
or a squashed [TU] spectrum. On the other hand, the stubborn agreement between SM 
predictions and observations in channels with large missing transverse energy (MET) cuts 
may indicate that the assumption of exact R-parity conservation is incorrect. 

Models with R-parity violation (RPV) have been proposed since the early days of 
SUSY pXJ. More recently a possible connection between the problems of baryon and lepton 
number violation and large flavor changing operators was highlighted [T2"j [T5]. The assump- 
tion of Minimal Flavor Violation (MFV) has been shown to be sufficient to prevent both 
rapid nuclear decay (and other baryon-number violating processes) and large corrections to 
flavor observables in the B, D, and K systems. In models of MFV SUSY, sparticles are 
pair produced as in R-parity conserving models, while the lightest supersymmetric partner 
(LSP) is generally unstable on collider scales and will decay via the baryon-number violating 
u c d c d c superpotential term. In this paper, we investigate one of the interesting scenarios 
that can arise in the model of [13]. 



In MFV SUSY models it is particularly compelling to consider the case when the LSP 
is a third generation squark: naturalness requires light third generation squarks in general, 
and we will see below that in the MFV scenario there is a high probability for this to be 
actually realized, due to the large top Yukawa coupling. The phenomenology of this scenario 
is also particularly rich, as the lifetime of an LSP stop or sbottom is long enough that the 
squark hadronizes to form a mesino by binding to a light quark pulled from the vacuum. It 
is, however, usually sufficiently short-lived to decay before reaching the detector. Observing 
squark production is challenging in this scenario, due to the lack of any obvious handles on 
the events, such as missing energy or displaced vertices. Instead we will make use of the idea 
of mesino-antimesino oscillations, following Sarid and Thomas [14J. We will demonstrate 
that sbottom-LSP pair production often allows for mesino-antimesino oscillations, which 
may lead to same-sign dilepton signals. 

A sbottom LSP decays dominantly to a top quark and a strange quark [13]. If one of the 
sbottom mesinos oscillates before decaying, the tops will be of the same charge, and if both 
tops decay leptonically this leads to same-sign leptons. These events would also contain b 
quarks from the top decays, providing further handles on the event. Recently, CMS searched 
for such events and placed bounds on their cross sections pj)]. We will show that this CMS 
search already places some bounds on sbottoms, which should improve significantly with 
more data. 

The rest of this paper is organized as follows. In Section 2, we study the typical squark 
spectra in MFV SUSY scenarios and demonstrate that the stop and sbottom are most often 
the lightest squarks. In Section 3, we present a calculation of the decay rate and oscillation 
time for a squark LSP in MFV SUSY and show that a significant oscillation probability is 
possible and occurs frequently. In Section 4, we comment on the sensitivity of existing LHC 
searches to this scenario. We conclude in Section 5. The details of the calculation of the 
mesino-antimesino oscillation rate is given in the appendix. 



II. MFV SQUARK SPECTRA 

In an MFV SUSY model the LSP decays and we are not restricted to models with a 
neutralino LSP, whereas the phenomenology of the model will depend on the identity of the 
LSP. In particular, the LSP can be colored and, as we consider below, can be a squark. 

MFV requires that all flavor violation be proportional to the appropriate combination of 
Yukawa matrices, which are treated as spurions of the flavor symmetry. The squark mass 
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FIG. 1. Distribution of lightest squark flavor over a random sampling of MFV SUSY parameter 
space. 



matrices are then required to have the following form [13 
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The D terms are automatically flavor diagonal and given by 
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where Q = +2/3 (—1/3) and T 3 = +1/2 (—1/2) for the up-type (down-type) squarks, s w is 
the sine of the Weinberg angle, tan is the ratio of the Higgs VEVs, and rriz is the mass 
of the Z. The parameters mf, cii, bi, and Ai arise from supersymmetry breaking, and we 
therefore expect them to be of order m 2 oft . 

Given these constraints, we can perform a scan over the parameter space that determines 
the squark spectrum. We select random values for the undetermined dimension-two parame- 
ters uniformly in [-m^, ^soft]- For this scan, we choose m so f t = 1 TeV and tan (3 = 10. The 
overall result is not very sensitive to this choice. We impose the constraint that the smallest 
eigenvalues of both squark mass matrices be greater than the top mass, mt ~ 175 GeV. In 
general, left-right mixing is not too large and we therefore use notation where bi refers to the 
mass eigenstate of the sbottom that is mostly a left-handed sbottom. We also impose that 
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the lightest stop-like squark be lighter than 500 GeV as demanded by naturalness. Under 
these conditions, the distribution of lightest squark flavors is given by Fig. [TJ 

We observe that roughly 85 % of parameter points have a third-generation lightest squark, 
out of which 15 % have a sbottom squark at the bottom of the spectrum. The large likelihood 
of a third generation lightest squark can be explained by the relatively large left-right mixing 
for this generation. This mixing tends to drive the mass of the lighter third generation squark 
down, making it more likely to be lightest overall. (There is also a significant contribution 
from the naturalness cut, since requiring one light stop tends to reduce the incidence of 
both stops being made heavy by a large positive a q .) Note that at large tan j3 this effect is 
enhanced for the sbottoms, making it even more likely to get a sbottom LSP. It is therefore 
natural to consider signatures of a sbottom LSP in MFV SUSY and we do so from this point 
on. 

III. MESINO OSCILLATION IN MFV SUSY 

The MFV SUSY scenario offers a rich phenomenology due to the naturally small decay 
width of the LSP, a consequence of approximate -R-parity conservation. The couplings are 
sufficiently small to yield LSP lifetimes that are longer than the timescales of SM short- 
distance physics, such as hadronization, yet often shorter than the timescales set by macro- 
scopic distances in the LHC detectors. In this intermediate range, it can be difficult to 
construct observables that are not overwhelmed by SM background. If the LSP carries 
color, however, then it lives sufficiently long to hadronize, an intriguing possibility. This 
process can yield additional phenomena that allow for efficient selection of SUSY events. 

The case of a sbottom LSP is particularly fruitful. If the gluino is heavy, the dominant 
SUSY production mode will be sbottom pair production. The dominant decay of the sbottom 
in MFV SUSY is to top and strange. The top has a leptonic decay mode, which already 
suppresses many SM backgrounds. As we show, the fact that the sbottom hardronizes allows 
for the possibility of sbottom oscillations, which lead, some fraction of the time, to same-sign 
lepton events. 

While other squark flavors can also oscillate, this turns out to be parametrically rarer. 
In addition, up-type squark LSPs do not decay leptonically, precluding the possibility of a 
same-sign dilepton signature. We do not consider these other possibilities further in this 
work. 

We also do not consider the case where gluino pair production is significant. This would 
lead to additional same-sign lepton events due to the Majorana nature of the gluino, provid- 
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FIG. 2. The leading diagram for the imparity violating sbottom decay. 



ing a background to the case we are considering. 

We begin by calculating the sbottom decay width. We denote the lightest sbottom mass 
eigenstate by b. Its decay width depends on an overall (generically order 1) coefficient that 
we denote by A". The Lagrangian terms that gives the decay is [13J: 

C = -(A ) tiju—Uq , — ±V fj — -huAd^ ~ -(A ) V td —^ M T s T , (4) 

where v = 174 GeV, V is the CKM matrix, and we use b\ to denote the lightest down-type 
squark, which we assume is predominantly sbottom-like. The mixing matrix U D is defined 
such that 

q q = u q %. (5) 

In this notation, q q are the squark flavor-basis fields in the mass basis of the quarks and qi 
are the squark mass-basis fields. The approximation is valid if the lightest sbottom is mostly 
right-handed. Otherwise, there is an additional suppression from the left-right mixing. The 
partial decay width can then be calculated using the diagram in Fig. [2j The result (neglecting 
the mass of the down quark in the phase space integral) is: 
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To gain some intuition about sbottom LSP decays, we now make a few approximations. 
The decay is dominated by b — > t c s c provided there is sufficient phase space. In the interesting 
segment of parameter space, the LSP is made up almost entirely of some admixture of the 
left-handed and right-handed sbottom, so that the decay width is approximately: 

2 

(2.6 x 1(T 10 GeV) | A" | 2 sin 2 # I ^ ) ( — ^— ) . (7) 
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FIG. 3. Diagrams for sbottom mesino oscillation mediated by a gluino. There are also similar 
diagrams mediated by neutralino exchange, which can become important if the gluino mass is very 
large. 

where 9 is the left-right mixing squark mixing angle. 

The sbottom decay rate is much less than the hadronization scale Aqcd ~ 0.2 GeV. 
Thus, the sbottom squark will hadronize before decaying to form fermionic mesino bound 
states B q = b*q and B c = bq c . If q = d, s, then the mesino is neutral, opening up the 
possibility for mesino oscillations, first discussed in [H] . Since few details of the calculation 
of the oscillation rate were given in I14j . we elaborate on it in Appendix |Al explaining the 



necessary approximations. Our final result, eq. (A17), is in broad agreement with that of [T 
and we restate it here: 
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This result depends on the nature of the spectator quark. We can use MFV to approx- 
imate the ratio of the oscillation rates as we have \Ul]\ oc |Vt g T4&| for M = L,R. In this 
approximation, we get: 
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The dependence of this ratio on the dimension-two parameters of the squark mass matrix is 
generically very weak. 

With this factor in mind, we consider oscillation of the sbottom-down mesino. The 
oscillation rate can be estimated by 
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~ (4 x 10- GeV) ( '* J ) cos= f 10 °° GeV ) . (10) 
v ; V 28 - 7MeV / V m 9 J 

These results are not too far from the decay rates, eqs. ([7]), but with different parametric 
dependence. Thus, we expect some parts of parameter space where the oscillation rate is 
comparable to or larger than the decay rate, leading to appreciable mesino oscillations. 
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FIG. 4. Oscillation parameter Xd = Am^/r and left-right mixing angle 9 resulting from a scan 
over parameter space, where 9 = corresponds to a pure left-handed LSP. 

To get a better sense of how common such a phenomenon is, we define the oscillation 
parameter 

x = —. (11) 

The time- integrated probability for a sbottom mesino to oscillate into an anti-sbottom mesino 
before decaying is 

**> s ^ = (12) 

The oscillation probability is small for x <C 1 and becomes appreciable near x ~ 1, whereas 
for i>1 the B oscillates very rapidly, and the mesino contains an equal mixture of sbottom 
and anti-sbottom components. We scan over parameter space using the same procedure as 
in Section |nj selecting points with a sbottom LSP and calculating Xd (the B& oscillation 
parameter) and 9 for each such point. The results of the scan are shown in Fig. |4j We 
observe that Xd > 1 in a significant portion of parameter space, particulary when the LSP 
is predominantly left-handed. 

If the sbottom is the LSP and has a mesino oscillation time comparable to or larger than 
its lifetime, then there is a very distinct signature of direct sbottom pair-production. The 
sbottoms will hadronize and the resulting mesino may oscillate before decaying. The mesino 
must be neutral for oscillations to occur, which occurs when the spectator is a down or 
strange quark, or roughly half the time as estimated from the B system. If exactly one of 
the mesinos oscillates before decaying, then the resulting two halves of the final state will 
have the same charge. Furthermore, these final states each involve a top quark whose charge 
is easy to tag if it undergoes a leptonic decay. This final state has same-sign leptons, b jets, 
and a small, but non-negligible, amount of missing energy. The entire chain is illustrated in 
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FIG. 5. Diagram for i?-parity violating sbottom decay that leads to same sign leptons. 
Fig. [5j The branching fraction for this mode is given by 

Br(bb* -> bbtt) = Br(W -> £u) 2 f(x d ,x s ) « f(x d ,x s ) x 6.5%, 

f(x d ,x s ) = 2 /kpfo) I 1 - M^i) I > ( 13 ) 

i=d,s \ j=d,s / 

where f(xd,x s ) denotes the probability that exactly one of the two mesinos oscillates. Here 
hi is the fraction of sbottoms that form mesinos with spectator i and we use the effective 
leptonic rate for the W which includes leptonic tau decays. Note that for Xi 3> 1 the rate is 
maximal, and since hd + h s w 1/2 we have f(xd,x s ) w 3/8. Despite the modest branching 
fraction, this decay mode will likely be the most sensitive channel for discovering a sbottom 
LSP in MFV SUSY. 
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FIG. 6. Recasted CMS bounds on sbottom direct production in terms of the sbottom mass and 
Xd- Only the four most sensitive signal regions are shown: SR0 in dashed green, SRI in dotted 
blue, SR2 in solid red, and SR4 in dash-dotted orange. The most conservative upper limit on the 
number of new physics events is used for each search region, though varying this number has little 
effect on the bounds. 

IV. BOUNDS FROM A CMS SEARCH 

CMS already has a search [15] that is quite sensitive to the above decay chain. The same- 
sign dilepton and b jets search includes search regions with 0, 30 GeV and 50 GeV MET 
cuts, all of which can be sensitive to our scenario due to the neutrinos from leptonic top 
decays. The relevant bounds from this search are presented in Table 2 of [15]. In addition, 
they present efficiency fits for the various cuts in terms of parton-level objects, allowing for 
easy reinterpretation. In this section, we use this information to reinterpret their bounds in 
terms of MFV SUSY with a sbottom LSP, and comment on future prospects. 

To obtain a bound, we generated pp — > bb* at 8 TeV using Pythia 8 with all showering 
and hadronization turned off. The events are decayed at the parton-level. The analysis cuts 
are applied using the efficiencies presented in Section 6 of [15]. No mixing is introduced 
in event generation, but an x^-dependent factor is applied to the final efficiency to account 
for the branching fraction to same-sign leptons. The cross-section for pair production is 
calculated at NLO using Prospino 2.1 JT6J [17]. The resulting bounds are shown in Fig. [6] 
SR2, which counts only positively chaged same sign pairs, yields the strongest bound, since 
a presumed fluctuation in the data led to all observed same-sign events having negatively 
charged leptons. For a maximal same-sign branching fraction, the exclusion extends between 
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180 GeV and 305 GeV. 

Note that obtaining same sign lepton events requires x > 1. The reason that there is 
sensitivity in the small x,i region is due to the possibility of producing strange mesinos. 
Even for x<i <C 1, it is possible to get x s > 1. 

V. CONCLUSIONS 

MFV SUSY is a compelling new paradigm for exploring supersymmetry without -R-parity 
that offers many new and challenging channels to explore at the LHC A systematic study 
of the phenomenology of all plausible scenarios in this framework is required to ensure full 
sensitivity to weak-scale supersymmetry. We have explored one interesting scenario with a 
sbottom-like squark LSP. 

Direct squark production will be essential for probing all possible corners of natural SUSY 
parameter space. Our work has demonstrated that the LHC can be sensitive to directly 
produced sbottom LSPs in the MFV SUSY scenario using the important fact that in this 
framework strongly-interacting LSPs will hadronize. Using a CMS search for same sign 
dileptons and 6-jets, we have put a bound on sbottoms with masses between 180 and 305 
GeV that undergo large mesino oscillations, which is a plausible scenario for the case with 
sbottom LSPs. 
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Appendix A: Determination of the Mesino Oscillation Frequency 

In this appendix, we present the details of the calculation of the mesino oscillation fre- 
quency, carefully listing all approximations as they enter. Throughout the appendix, we will 
assume a down quark spectator, but the results extend trivially to the strange quark case. 
We further denote the lightest squark by b and assume that it is sbottom-like. 

In the quark and squark mass basis, there are two combinations of sbottom and down 
quark that correspond to light mesino Weyl fermions: 

B l = B = b*d, B 2 = B C = ~bd c . (Al) 
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The most general quadratic Lagrangian for these mesino fields is given by: 

£ = ^m ij B i B j + h.c. . (A2) 

Before including corrections due to the gluino, the diagonal entries of m^ vanish, and the two 
Weyl fermions combine to form a Dirac fermion. The mass, corresponding to the off-diagonal 
terms in ( A2 ) , is given to leading order by 

= m~ b ■ (A3) 

The leading corrections are of order Aqcd, which we neglect. 

The diagonal elements mn and m 22 , corresponding to Majorana masses for B and B c , are 
not in general equal, and are generated at leading order by tree- level gluino exchange, leading 
to an oscillation between mesinos and antimesinos. The oscillation frequency is equal to the 
mass splitting between the two mass eigenstates, whose squared masses are the eigenvalues 
of rn^rn. We take myi to be real by performing an appropriate field redefinition, in which 
case the eigenvalues of m)m are given by 

^ ^|mn| 2 + \m 22 \ 2 + 2mj 2 ± \J (|mn| 2 + \m 22 \ 2 + 2m\ 2 ) 2 ~ M m l 2 ~ m n m 22\ 2 ^j ■ (A4) 
To leading order in mu and 77722, the resulting mass splitting is 

uj = Am = \m u + m^l • (A5) 

We work at leading order in the heavy squark approximation. Instead of determining mu 
and 777 2 2 directly, we employ the simple and general formula: 

u = —\(B(d,s)\U cS (d)\B(d,s))\, (A6) 

777 x2 

for uj 777 12 , where W e ft(x) is the effective Hamiltonian density generated by integrating 
out the gluino and \B(p,s)) and \B(p, s)) denote one-particle mesino and antimesino states, 
respectively, with momentum p and spin s with no sum over s. (We use the standard 
covariant normalization for one-particle momentum eigenstates, (p\q) = 2E^(2tt) 3 S^ (p — q) .) 
The effective Hamiltonian density from integrating out the gluino is: 

Heff = ^(b*d)(b*d) + ^(bd c )(bd c ) + h.c. , (A7) 

for coefficients Cl and Cr to be determined, where the color indices are contracted as 
indicated by the parentheses. Thus, 
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^(B\{b*d){b*d)\B) + ^(B\(b*d c ^)(b*d^)\B) 



(A8) 
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The structure is very similar to (A5), and indeed the two terms within the absolute value 



in (A8) are precisely times the Majorana masses which appear in (A5). 

To determine the Cl r, we compare the short-distance amplitudes for oscillation obtained 



using the MSSM Lagrangian and using the effective Hamiltonian in (A7). The MSSM gluino 
exchange amplitudes Ml and Mr (Fig. |3]) are given by 
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where we work in a basis where the gluino mass is real, and an overall factor of two arises 
since the gluino-quark-squark vertex comes with a factor of \/2. The color factors in these 
amplitudes simplify to [18J: 



j.a j.a _j_ j.a +a 
ij i'j' ' ij' i'j 



The effective operators in (A7) yield amplitudes: 

•M-L R = CL,R,(8ij>8i'j + SijSi'j') . 



(A10) 
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By demanding that Ml (Mr) from (A9) is equal to M' L (M' R ) from (All), we extract the 
coefficients Cl and Cr\ 
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(A12) 



The same result can be obtained in the large rrig limit by integrating out the gluino in the 
Lagrangian, neglecting the kinetic term. 



As QCD is parity invariant, the hadronic matrix elements in (A8) are equal. We esti 



mate them using the vacuum insertion approximation. In this approximation, we insert the 
vacuum between the operators in all possible ways, giving [191 120] : 



{B\(b*tf){b*dt)\B) « 2 {B\{b*^)(Q\{bW)\B) + (Bmd>)\0)(0md*)\B) 



(A13) 



where we indicate color indices explicitly, and there are two ways to obtain each of the 
terms, yielding a prefactor of 2. The contraction with the color-neutral external state kills 
the terms with % ^ j in the second term. Exactly one in every A^ c terms has % = j, so we get 
the relation: 

(B\(b*dn\0)(0\(b*d*)\B) = ^(M*d^\0)(0\(b*d^\B). (A14) 
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Our result is thus: 



(B\(b*d)(b*d)\B) = {B\(b*d c] )(b*d c] )\B) 

!^^(S|(6*d*)|0>(0|(6^)|B> =2^^/|m 5 . (A15) 



The mesino decay constant can be estimated using the B meson decay constant and 
assuming heavy quark symmetry. Up to threshold corrections, the relationship is given by 
J2TJ: 

= rw f ^MV /23 ( 2^L\ m . ( Ai6) 

V m ' bl \a s {jn t )J \a s {m- bi )J 



Using the latest values of f B = 190.6 MeV [221 123] and a s (m z ) = 0.1184 [23], the MS quark 
masses m;, = 4.19 GeV and mt = 160 GeV [23], as well as with a numerical solution to the 



NNNLO beta function for a s [25], which we evaluate at mg = 300 GeV, we find a value of 
/_§ = 28.7 MeV. 

Putting these pieces together, we arrive at our final expression: 

A- = 9l IK,) 2 + «i) 2 | /IK) (l " ^) , (A17) 



with fg given by (A16). This agrees with [T3] up to a factor of 8 and the dependence 
on the CP-violating phase in the squark mixing matrix. This result has some hadronic 
uncertainty, which we estimate to be of order 10 % based on estimates of the validity of the 
same approximations for the B meson systems. 
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